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ABSTRACT 

SONYC - Suhstellar Objects in Nearby Young Clusters - is a survey program to investigate the 
frequency and properties of suhstellar ohjects with masses down to a few times that of Jupiter in 
nearby star-forming regions. In this third paper, we present our recent results in the Chamaeleon- 
I star forming region. We have carried out deep optical and near-infrared imaging in four bands 
{I,z,J,Ks) using VIMOS on the ESO Very Large Telescope and SOFI on the New Technology 
Telescope, and combined our data with mid-infrared data from the Spitzer Space Telescope. The 
survey covers ~0.25deg^ on the sky, and reaches completeness limits of 23.0 in the I-band, 18.3 in 
the J-band, and 16.7 in i^g-band. Follow-up spectroscopy of the candidates selected from the optical 
photometry (/ < 21) was carried out using the multi-object spectrograph VIMOS on the VLT. We 
identify 13 objects consistent with M spectral types, 11 of which are previously known M-dwarfs with 
confirmed membership in the cluster. The 2 newly reported objects have effective temperatures above 
the suhstellar limit. We also present two new candidate members of Chamaeleon-I, selected from our 
JK photometry combined with the Spitzer data. Based on the results of our survey, we estimate that 
the number of missing very-low-mass members down to ~ O.OOSMq and Ay < 5 is < 7, i.e. < 3% of 
the total cluster population according to the current census. 

Subject headings: stars: circumstellar matter, formation, low-mass, brown dwarfs - planetary systems 



1. INTRODUCTION 

Brown dwarfs {M < 0.08 Mq) are an important class 
of object to test the mass dependence in the formation 
and early evolution of stars and planets. As of today, 



the origin of these suhst ellar objects it not clear (Whit- 
worth fc Goodwin| | 200 5 [) ; they m ay be the result ot tur- 



bulent fragmentation ( [Padoan fc Nordlund 2004]) or ejec 
tions from multiple systems ( [Bate||2009| ). In addition, 
brown dwa rfs might have been eiectedlrom protoplane- 
tary disks ( Stamatellos fc Whitworth]|2009| ) where they 
have formed by disk fragmentation, and we might ob- 
serve the overlap of star and planet forming processes. 
Further issues that can be addressed by investigating 
young brown dwarfs are, among others, the evolution 
of accretion disks (e.g., the impact of photocvaporation) , 
the formation of multiple system s (jB atc 2009 ) , the mass 
dependence of planet formation ( [Pay ne & Lodato 2007), 
and the ea rly regulation of angular momentum (iScholzj 
fc Eisloffel 2004, 2005). 

'I'he prerequisite for tackling these issues is to ob- 
tain a reliable picture of the frequency and properties 
of young brown dwarfs. The currently available sur- 
veys are affected by two major biases: a) In most re- 
gions they do not extend below the deuterium-burning 
limit 0.015 M0), i.e. they do not cover the regime of 
free-floating objects with planetary masses. In several re- 
gions, however, it ha s been shown that the mass function 
extends to 0.01 (Zapatero Osorio et al. 2000 Lucas 
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et al. 2005 Luhman et al. 2008) or even below, and a 
lower mass limit (the "bottom of the IMF" ) has not been 
found yet, which indicates that our census of suhstellar 
objects is incomplete at the lowest masses; b) Some of the 
deepest surveys are biased by design, either because they 
are based on mid-infrared data from Spitzer and th us will 
only find objects surrounded by dusty disks (e.g., Luh- 



man et al. 20081 or because they are carried out usmg 
narrow- band methane imaging, whi ch probes only a spe 



cific range of spectral types (e.g., Haisch et al. 2010). 
While samples obtained with these methods are useful 
to constrain the mass function, they are of limited use 
for other purposes (for example, Spitzer samples cannot 
reliably probe disk lifetimes as a function of mass) . 

Our program SONYC (Suhstellar Objects in Nearby 
Young Clusters) aims to mitigate these selection effects 
and to help in establishing a more complete census of 
brown dwarfs. The project is based on three basic prin- 
ciples: 

1) By making extensive use of 4- to 8-m-class telescopes, 
we aim to reach mass limits of 0.005 Mq and below, with 
the specific goal of probing the bottom of the IMF. 

2) The primary means of identifying candidates is broad- 
band imaging in the optical and near-infrared, thus aim- 
ing to detect the photosphere. This results in large sam- 
ples of candidates, i.e. requires extensive follow-up spec- 
troscopy. This dataset is complemented by publicly avail- 
able Spitzer data. 

3) By probing several star forming regions we want to 
probe for environmental differences in the frequency and 
properties of suhstellar objects. 

So far, we have published results for two r egions, 
NGC 1333 jScholz et al.|[2009l ) and p Ophiuchus (iGeers 
et al.||2010j) In NGC1333 we identify 19 substellaFoE 



jects, 12 of them previously unknown. We tentatively 
find a cutoff in the mass function between 0.01 and 
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Fig. 1. — Projected distribution of the observed fields in Cha-I. Solid lines show four fields observed by VIMOS in the imaging mode. 
Each field consists of four quadrants; labeling of the VIMOS fields (Fl to F4) is shown near the centers of the gaps between the quadrants. 
Spatial distribution of the SOFI fields approximately matches the VIMOS fields. IRAC fields are outlined with dashed (deep) and dash- 
dotted line (shallow). Optical spectra wer e obtained with V IMOS in the MPS mode, for the objects in fields Fl, F3, and F4. Black dots 
mark the known members (Luhman 2007, |Luhman fc Muench ^ 2008 .Luhman et al.|2008| l, with brown dwarfs (T^ff < 3000 K) additionally 
shown as diamonds. Squares mark tke two candidates trom iablej^ and red crosses mark the two new candidate members from Table [2] 

0.02 Mq. This may be the first sign that we have reached 
the bottom of the IMF. In p Ophiuchus we spectroscop- 
ically confirm a previously unknown brown dwarf, and 
identify a number of new candidates based on Spitzer 
data. In this region, however, the current census is highly 
incomplete due to strong and variable extinction. 

In this paper we present the SONYC results for the 
star forming region Chamaeleon I (hereafter Ch a-I). Cha- 



I har bors a rich population of known YSOs (Luhman 



20071, most of them grouped in two clusters, each dom- 



inated by a bright A star. Cha-I is less compact than 
NGC1333 and has less extinction than p-Oph. Pre- 
vious surveys hav e identified a significant number of 
substellar obje cts ( jComeron et a l. 1999 Neuhauser fc 
Comcron 1999; Oasa et al. 1999; C Weron et al.||!^004 ' 
Eopcz Marti et al.,|2004J . Extensive photometric anc 



spectroscopic studies of the region jLuhman et al.||2004 



Luhman|[2007| [Luhman fc MuenchPOOS, |Luhman et al 



2008p result in 40 spectroscopically contirined objects be- 
low the hydrogen-burning limit. Among these are a few 
objects with es timated masses at or below the deuterium- 
burning hmit dLuhman et~aI1 [2b06 , 2008). Throughout 
this paper we will use the current membership census, 
referred to as the "Luhman census", by Luhman (2007, 
215 ob .iects, his table 6), co mplemented by the n ew mem- 
bers in Luhman fc Muench (2008, 8 objects) and Luhman 



et aL| ( |2008[ 4 objects). 

'i'his paper is organized as follows. Observations and 
data reduction are described in Section [21 In Section |3] 
we describe selection of substellar candidates from pho- 
tometry, and the spectroscopic follow-up of the candidate 
objects in Cha-I. The results are discussed in Section [Ij 
and conclusions are presented in Section [Sj 

2. OBSERVATIONS AND DATA REDUCTION 

This study is primarily based on data from the Euro- 
pean Southern Observatory, obtained in runs 078.C-0049 
and 382.C-0174. The first campaign provided deep opti- 
cal and near-infrared imaging. From the optical dataset 
we selected a list of candidate members which were ob- 
served using multi-object spectroscopy in the second run. 
Additionally, we made use of the archival Spitzer images 
for our target regions. In the following we describe the 
observations, the basic reduction and photometry. The 
fields covered by our surveys are outlined in Figure [ij 

2.1. Optical imaging 

We used VIMOS at the VLT to observe Cha-I in the 
optical wavelength regime. VIMOS is a wide-field imager 
with 4 CCDs arranged in a 2 x 2 array. Each detector 
covers a field of view of 7' x 8' with a pixel resolution of 
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Fig. 2. — Histograms of the objects in our (a) optical and (b) near-infrared photometric catalogs. The peaks at 7 = 23.0, J = 18.3, and 
Ks = 16.7 indicate the completeness limit in each band. The faintest objects in the survey are found at 7 > 25, J > 19.5, and Kg > 17.5. 



0'.'2. The four quadrants are separated by gaps of approx- 
imately 2'. Four pointings were observed, each in the I- 
and z-band filters. See Table [l] for more details about 
the VIMOS imaging run. With this observing strategy 
we do not cover the gaps between the chips, i.e. the sur- 
vey is spatially incomplete. Most of the currently known 
young population of Cha-I is clustered around the two 
bright stars HD 97048 and HD 97300. These two areas 
(with the bright stars in the central gaps of the VIMOS 
array) are covered in our fields 1 and 4. The two remain- 
ing fields are located adjacent to each other and north of 
field 1. Taken together, the SONYC fields cover 38% of 
the known Cha-I population, including 17 known brown 
dwarfs. 

A standard image reduction was preformed for each in- 
dividual science observation, using the VIMOS pipeline 
recipes. The individual images were then combined and 
the world coordinate system was calibrated against ob- 
jects found in the 2MASS point source catalog with 
17 > J > 12. 

Sources were identified and extracted from each CCD 
chip usin g the Source Extractor {SExtractor) software 
package (Bertin & Arnouts 1996). For object identifi- 
cation we required at least 5 pixel to be above the 3cr 
threshold detection limit. SExtractor automatic aper- 
ture fitting photometry routines were used to calculate 
the flux of each detected source. We further rejected 
objects with flux signal-to-noise of less than 5, overly 
elongated objects [a/h < 2.0), and objects which could 
not be found in both the / and z catalogs. 

To calibrate the I-band photometry we compared our 
instr u menta l magnitudes with those in Table 1 of |Luh- 



man 



(2007), for objects with / > 16. No photometric 
calibration was available for z-band, so in order to com- 
pare color magnitude diagrams for different CCD chips 
and fields, we calculated the median instrumental I — z 
of all catalog objects in the range 20.0 < I < 17.5, 
— 1.5 < I — z < 1.0. We then subtracted the difference 
between the median I — z for each chip and the median 
for F2, CCD chip #1. 
From the histogram of the magnitudes (see Figure [2]), 



we derive the completeness limit of the optical survey at 
/ — 23.0 ±0.2 mag, and the limiting magnitude of about 
25 mag. 

2.2. Near-infrared imaging 

The near-infrared observations were designed to pro- 
vide J- and i^s-band photometry for the same area cov- 
ered in the optical survey with VIMOS. We used SOFI at 
the NTT in its wide-fleld mode resulting in a field of view 
of 4.9' X 4.9' and a pixel scale of 0'.'288. The definitions 
o f the SOFI filte rs are comparable to the 2MASS bands 
(|Carpenter|2001[ ). See Table [l] for more details about the 
S01<'1 run. 

Standard data reduction including sky subtraction, 
bad pixel, fiat field, and the detector row cross-talk cor- 
rection, was performed using IDL and DPUSER[^ (T. 
Ott). The 32 integrations per VIMOS quadrant were 
combined into a single mosaic using a simple shift-and- 
add algorithm. Shifts between the individual exposures 
were determined by cross-correlation, performed by the 
jitter rou tine (part of the ESO eclipse package; Devil- 
lard 1997). The images were then shifted and rnedian 
averaged. 

For the source detection we used SExtractor, requiring 
at least 5 pixels with fiux above the 3.5 a detection limit. 
Saturated and extended sources were rejected from fur- 
ther consideration. For the photometric calibration, we 
identified all 2MASS sources matching the J- and Ks- 
band sources in our images with signal-to-noise better 
than 10. The number of useful calibrators per field ranges 
from 85 to 164. With these stars we determined the ze- 
ropoint ZPl and ZP2 and a color coefficients CI and C2 
to transfer our magnitudes (j, k) to the 2MASS system: 
J = j + ZPl + Cl{j -k) &nAK = k + ZP2 + C2{j - k). 
With one exception, the color terms are 0.02 or smaller, 
and do not have significant effects on the calibration, 
which is in agreement with the flndings in Carpenter 
(2001). The zeropoints are 1.15 to 1.23 mag in J and 
-0.53 to -0.59 mag in Ks- The spread can be explained 
with the variations in airmass between our fields. The 

^ http:/ /www. mpe.mpg.de/~ott/dpuser 
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TABLE 1 

Observing log for our imaging and spectroscopic surveys. 





(Xyo jii\J\j\j J 


x( T9^^^^ 

J ZUUU J 




UT date 


airmass seeing^ 


notes 










Optical and near- 


infrared imaj 


;ing 


1 


11:07:30 


-77:39:00 


I, z 


Jan 17, 2007 


1.7 


0.9" 


11 X 150 s in z, 10 X 150 s in / 


2 


11:05:00 


-77:21:00 


I, z 


Feb 10, 2007 


1.7 


0.7" 


11 X 150 s in z, 10 X 150 s in I 


3 


11:10:00 


-77:21:00 


I, z 


Feb 11, 2007 


1.7 


0.8" 


10 (9 for Q3)xl50s in z, 11 (10 for Q4) xl50s in / 


4 


11:11:00 


-76:36:00 


I, z 


Feb 05, 2007 


1.6 


0.5" 


10 X 150 s in z, 10 X 150 s in I 


1 


11:07:30 


-77:39:00 


J, Ks 


Nov 11, 2006 


1.9 


1.0" 


64 X 20 s in J, 64 X 8 s in Ks 


2 


11:05:00 


-77:21:00 


J, Ks 


Dec 05, 2006 


2.0 


1.2" 


64 X 20 s in J, 64 X 8 s in Ks 


3 


11:10:00 


-77:21:00 


J, Ks 


Dec 20, 2006 


1.6 


0.9" 


64 X 20 s in J, 64 X 8 s in Ks 


4 


11:11:00 


-76:36:00 


J, Ks 


Nov 16, 2006 


1.9 


1.1" 


64 X 20s in J, 64 X 8s in Ks 












Spectroscopy 




1 


11:07:30 


-77:39:00 


LR_red 


Feb 16, 2009 


1.7 




4 X 1610s 


3 


11:10:00 


-77:21:00 


LR_red 


Mar 05, 2009 


1.7 




4 X 1610s 


4 


11:11:00 


-76:36:00 


LR_red 


Mar 06, 2009 


1.7 




4 X 1610s 



^ measured as a typical FWHM of unsaturated sources in our images 

typical error in the photometric cahbration is 0.06 mag. 
From the histogram of the magnitudes (see Figure [2]), we 
derive the completeness limits of the NIR part of the sur- 
vey: J = 18.3 ±0.2 mag, and Ks = 16.7 ±0.2 mag. The 
limiting magnitudes are about J = 19 and Ks = 17.5. 

The WCS system for all images was determined with 
ccmap in IRAF, based on the 2MASS sources. Typically, 
we used 50-150 2MASS stars per VIMOS quadrant; the 
resulting accuracy in RA and DEC is ^ 0'.'15. 

2.3. Spitzer photometry 

In addition to our optical and near-infrared images, we 
use IRAC data from the Spitzer Space Telescope to iden- 
tify candidates. The regions of our optical/near- infrared 
survey are covered by five AORs: 3955968, 6526208, 
12620032 from the program #36 (PI: Fazio), 3651328 
from program #6 (PI: G. Fazio), and 3960320 from pro- 
gram #37 (PI: G. Fazio). The three first mentioned 
AORs are much deeper than the other two (100 s vs. 12 s 
integration time). Our analysis is based on the pipeline 
reduced PBCD images in the four IRAC channels at 
3.6, 4.5, 5.8, and 8.0 /im from these five AORs. Us- 
ing daof ind within IRAF we searched for objects with 
peak fluxes exceeding 5 times the noise in the image. 
This yields signiflcant contamination, particularly in the 
bright cloud cores of Cha-I. Therefore we required an ob- 
ject to be detected in all four IRAC images of the same 
AOR within ±1 pixel. Note that the PBCD images for 
these AORs are sampled to 0.6" /pix. For the resulting 
clean object catalogs we carried out aperture photome- 
try using daophot. We chose relatively small apertures 
of 5 pix and a background annulus of 5-10 pix. We de- 
termined aperture corrections using a small sample of 
isolated bright stars, to shift our photometry to the stan- 
dard aperture of 10 native pixels, as given in the IRAC 
handbook. The IRAC photometry catalogs were cross- 
correlated with the object catalog from our deep J- and 
Ks-hand images. In this combined near/mid-infrared 
catalogs we will search for objects with the typical mid- 
infrared excess as expected for sources with disks. 

2.4. Multi-object spectroscopy 



The optical spectra have been obtained using VI- 
MOS/VLT in the Muhi-Object Spectroscopy (MOS) 
mode. The FOV in the MOS mode is the same as in 
the imaging mode. We have observed three out of four 
fields covered by the imaging run (Fl, F3, and F4), using 
the low resolution red grism (LR_red). We covered the 
wavelength regime between 5500 A and 9500 A, with the 
spectral resolution of R~210, slit width of 1" and a spa- 
tial scale of 0'.'205. The total on-source time for each field 
was 4x 1610s. Data reduction was performed using the 
VIMOS pipeline provided by ESO. Data reduction steps 
include bias subtraction, flat-field and bad-pixel correc- 
tion, wavelength and flux calibration, as well as the final 
extraction of the spectra. As no jittering offsets have 
been applied between the individual exposures, fringe 
correction could not be performed. Our spectra there- 
fore suffer strong fringing effects longwards of ~ 8000 A. 
In total, we have arranged 308 slits across the three ob- 
served fields (126 in Fl, 81 in F3, and 101 in F4). After 
discarding empty slits and very noisy spectra, we are left 
with spectra of 229 objects ready for inspection. 

3. SELECTION OF SUBSTELLAR CANDIDATES IN CHA-I 

3.1. Optical color-magnitude diagram 

The (I, I-z) color-magnitude diagram of the sources 
observed in the four Cha-I fields is shown in Figure [3| 
Previously spectroscopically confirmed M-type members 
(Luhman census) are shown as diamonds. On the right- 
hand side of the diagram, we plot the expected /-band 
brightness of Cha-I members with the mass 0.075, 0.015, 
and 0.005 Mq, each at Av= and 5. For the calcula- 



tion we use COND03 jBaraffe et al.|2003[ ) and DUSTYOO 
( Chabrier et al.||2000 ( evolutionary models. 

Very- low-mass (VLM) objects are expected to occupy 
a distinct region on the red side of the broad cumulation 
of the background main-sequence stars. We construct 
the selection box in the CMD such that it encompasses 
the majority of the known members of the spectral type 
M6 or later. 142 sources found within the selection box 
served as candidates for spectroscopic follow-up with VI- 
MOS. 

The faintest confirmed member in the color-magnitude 
diagram has / = 24.6. We identify this source as 
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Fig. 3. — I vs I — z color-magnitude diagram of the four Cha-I fields observed with VIMOS. The plus signs mark the substellar candidate 
sample, and small blue squares the objects with available optical spectra. Two M-dwarf candidates identified in this work are shown as red 
circles. Diamonds mark the confirmed M-type members objects according to the Luhman census. On the right-hand side of the diagram, 
we plot the expected /-band brightness of Cha-I members with the mass 0.075, 0.015, and 0.005 Mq, each at Av=0 and 5; based on 
COND03 and DUSTYOO evolutionary models. 



2M J11070369-7724307. [Luhman et al | ( |2008j^ published 
the spectrum, and estimated spectral type M7.5 ± 1 and 
effective temperature of 2795 K. The source s uffers the 
extinction of Ay ~ 15 — 16 (estimate from 'LuhmanI 
et al. 2008 and our J ~ K color, assuming the intrin- 
sic J — K = 1), which is higer than the average extinc- 
tion in Cha-I (see Section lij]) . According to DUSTYOO 
and COND03 models, a 2330 K (-M8.5) object at the 
distance of Cha-I and Ay — 16 is expected to have 
/ = 23.7 — 24.8, in agreement with our photometry. 

3.2. Spectroscopic follow-up 

Slits were allocated to 60 out of 142 objects selected 
from the optical CMD, using the VIMOS Mask Prepa- 
ration Software (VMMPS). A number of slits were also 
placed on the already known M dwarfs in the region. To 
fully use the multiplexing capability of the spectrograph, 
remaining slits were placed on the randomly chosen ob- 
jects within the FOV. As already mentioned above, some 
slits appear to be empty, or have poor signal-to-noise ra- 
tio. Blue squares in Figure [3] show 229 objects with the 
VIMOS spectra that can be used for spectral classifica- 
tion; 18 of those are found within the selection box. The 
approximate cutoff for spectroscopy is evident from the 
CMD, as the spectra of the vast majority of the objects 
fainter than / = 21 exhibit S/N too low to allow conclu- 
sive classification. 

Our goal is to identify young members of Cha-I with 
masses in the substellar regime, which typically cor- 



responds to objects with effective temperature below 
3000 K. The photometrically chosen sample might be 
contaminated by embedded stellar members of Cha-I, 
reddened background M-type stars, and (less likely) late 
M- and early L-type objects in the foreground. The red 
portion of M-dwarf spectra is dominated by molecular 



features (mainly TiO and VO; |Kirkpatrick et al. 1991 



1995 1, which allows a relatively simple preliminary selec 



tion of candidates based on visual inspection. Here we 
included all the objects with a clear evidence, but also, 
for the sake of caution, those objects with only a ten- 
tative evidence of molecular features. Another property 
that was used in the selection is the slope, which should 
be positive in the red portion of the optical regime for 
M-type stars. Earlier stellar types typically have neg- 
ative slopes, or appear flat (late K-type). Among the 
rejected spectra we mainly find flat spectra with no, or 
with very weak indication of molecular features (late K- 
dwarfs, or reddened earlier- type stars), and reddened fea- 
tureless spectra (background or embedded stars of spec- 
tral type earlier than K). 

After this preliminary selection, our sample contains 
28 spectra. To determine the effective temperature of 
our candidate objects we per formed spectral fit ting us- 
ing the AMES-Dusty models ( |Allard et al.|2001[ ). Model 
spectra are previously smoothed by a boxcar function to 
match the wavelength step of the VIMOS spectra. We 
compared the effect of log 5 in range 3.5-6.0 for constant 
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Fig. 4. — Spectra of the objec ts in Table p| ("black ) , and the best- 
fit model (blue; AMES-Dusty; |Allard et ar]|2001[ |. The spectral 
resolution of the models has been reduced by boxcar smoothing, in 
order to match the resolution of our spectra. Each spectrum has 
been corrected for extinction, but not for the atmospheric absorp- 
tion. 



effective temperature and found that gravity does not sig- 
nificantly change the shape of the optical low-resolution 
spectra. We search for the lowest by varying the ef- 
fective temperature Tes, and the interstellar extinction 
Ay, for the lowest available logg (3.5 or 4.0). We vary 
Teff between 2000 K and 3900 K, in steps of 100 K, and 
Av from to 6, in steps of 0.5 mag. To de-redden our 



spectr a, we apply the extinction law from Cardelli et al 
( 1989 1, assuming Ry = 4. Due to the low resolution, it is 



the overall shape and the broad features in the spectrum 
that finally determine the best fit. To keep the fitting 
range as broad as possible, we set the wavelength lim- 
its at 6200 and 8500 A. The cut at the red side is set 
to exclude the fringe-dominated portion of the spectra, 
while the cut at 6200 A is set to avoid the problematic 
region around 6100 A, where the models typically over- 
estimate the observed flux. The results of our spectral 
fitting are given in Table [2] and shown in Figure [4j We 
show only the objects with the best-fit Tgff at or oelow 
3800 K because the objects found at the edge of the fit- 
ting grid could be even warmer. The uncertainty is 100 K 
for the derived effective temperature, and 0.5 mag for the 
extinction, which reflects the spacing of the grid. 

The models predict a much sharper peak at 
^7000 A than apparent in our spectra, a discrepancy 
that become s more pronounced for l ower temperatures. 
According to Kirkpatrick et al. ( 1991 1, the TiO bandhead 
at 7050 A is found in all late-K to lat e-M spectra , but is 



strongest in early- to mid-M objects. Reid et al. ( 1995 1 
suggest that the depth of this feature ]Drovides a reliable 
spectral typing scheme for K7 - M6.5 dwarfs. At later 
types, however, the depth of the TiO feature decreases, 
which reflects both saturation in the TiO bandheads , and 



the presence of increasingly strong VO bandheads ( Reid 
et al. 1995). This behaviour is observed in our spec- 
tra, but is not well reproduced by the models. In the 
same spectral region, object #4 shows flux levels above 
the model prediction and in this aspect stands out from 
the rest of the spectra in Figure[4] This might indicate 
that the object is hotter than our oest-fit temperature at 
3800 K. 

We have to point out that the spectra shown in Fig- 
ure |4] are not corrected for telluric absorption. However, 
the two above mentioned spectral regions where models 
do not fit the spectra well (6100 and 7050 A) are free from 
atmospheric bands. The only region where the effects of 
telluric absorption are evident, is a deep O2 absorption 
at - 7600 A. 

The majority of the sources (11/13) shown in Figure |4] 
are previously identified members of Cha-I according to 
the Luhman census. The effective temperatures obtained 
by our model fits are in agreement with those from the 
literature. The mean of the difference between the two 
sets of Toff is 20 K, with a standard deviation of 150 K. 
The two newly discovered objects with M-dwarf spectral 
features have effective temperatures of 3400 and 3800 K, 
placing them among the early M-type objects. The ex- 
tinction estimate for these two sources lies in the range 
Ay = 1 — 3 which is consistent with the extinction to- 
wards Chamaeleon. However, without other clear signa- 
tures of youth, we cannot yet claim cluster membership. 
Further spectroscopy is needed to clarify their status. 
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TABLE 2 

Parameters of Cha-I M-type objects observed with VIMOS 



ID 


a(J2000) 


<5(J2000) 


Tcft (K) 


Av (mag) 


Reference^ 


name 


1 


11 09 45.1833 


-77 40 33.4442 


3000 


0.5 


L07 


tbO H« 566 


3 


11 08 16.9700 


-77 44 11.6052 


3000 


2.0 


L04 


ChaHal3 


4 


11 10 18.7603 


-77 46 13.8007 


3800 


4.0 


this work 




6 


11 10 11.5437 


-77 33 51.9489 


3400 


1.5 


L07 




10 


11 12 02.7459 


-77 22 48.1824 


2700 


0.5 


L07 




11 


11 10 50.7099 


-77 18 03.2794 


3300 


2.0 


L07 


ESOHa568 


12 


11 07 59.9412 


-77 15 31.9702 


3000 


3.0 


L07 


ESO Ha 561 


13 


11 08 55.9646 


-77 27 13.2001 


3000 


5.5 


L07 


ISO 167 


23 


11 13 01.2000 


-76 32 54.0509 


3400 


1.0 


this work 




25 


11 08 51.7863 


-76 32 50.4529 


2700 


2.0 


L04, L07 




26 


11 08 49.5341 


-76 38 44.0479 


2600 


4.0 


L07 




27 


11 09 52.1734 


-76 39 12.6947 


2800 


3.0 


L04 


ISO 217 


28 


11 10 06.5998 


-76 42 48.4937 


2500 


2.0 


L07 





'igurc[4| 

References for the published spectra: L04: [Luhman et al.| ( |2004[ |, L07: |Luhman| l |2007[ 




2 3 
J— K (mag) 

Fig. 5. — Near-infrared color-magnitude diagram for the sources 
in Chal. Known members are marked with crosses; objects with 
Spitzer counterpart are shown in red. The solid line shows the 
approximate limit of the combined JK — Spitzer catalog. 



X X 

X 

\>o< 



0.5 1 1.5 

13 - 14 (mag) 

Fig. 6. — IRAC color-magnitude diagram for the sources detected 
in J- and A'g-band. Crosses show 20 objects with color excess typ- 
ical for objects with disks and/or envelopes. Two of these sources 
do not yet have clarified membership (squares). 

3.3. Candidates from near-infrared and Spitzer 
photometry 

The second candidate selection was carried out based 
on our near-infrared catalog, which was c omp lemented 



by Spitzer/IRAC data, as outlined in Sect. 2.3 In total 



64 of them are known members of Cha-I, according to the 
Luhman census. 

In this catalog are 643 objects with Spitzer/IRAC 
counterpart within 3'.'0, most of them (521) in the deep 
images from program #36. In Figure [5] we show the 
(J, J — K) color-magnitude diagram generated from our 
catalog. Known Cha-I members are marked with crosses; 
objects with Spitzer counterpart are plotted in red. As 
can be seen in the diagram, the Spitzer survey misses 
many objects at the faint end of our JK catalog. The ap- 
proximate magnitude limit of the combined Spitzer— J K 
catalog is shown with the solid line. The limit ranges 
from J — 17 — 20 for Ay = 0—10 mag. According to 
the COND03 and DUSTYOO isochrones, this limit corre- 
sponds to masses of ~ 0.005 Mq. 

To select new candidate members in Cha-I, we make 
use of the standard IRAC color-color diagram, as shown 
in Figure [6| Plotted are the 5.8 — 8.0 /im color vs. the 
3.6 — 4.5 /im color for all 643 objects with Spitzer coun- 
terpart. As shown extensively in the literature, this di- 
agram serves as a robust diagnostic for the presence of 
circumstellar dust around a p articular object (|Allen et al 



2004| Hartmann et al. 2005). Objects with disks and/ or 



envelopes show excess m the two IRAC colors. The over- 
whelming majority of our sources appears around the ori- 
gin of the plot, which precludes the presence of a disk. A 
small group of 20 objects has significant excess in both 
axes and is clearly distinct from the bulk of the dat- 
apoints (marked with crosses). 18 of them are known 
members according to the Luhman census. The remain- 
ing 2 are listed in Table [3j and marked with squares in 
Fig. |6] To our knowledge, these objects are new candi- 
date members of Cha-I, although they both look slightly 
extended in the near-infrared images and could be galax- 
ies. For a definite statement about the membership of the 
new candidates, spectroscopy is needed. 



4. DISCUSSION 



4.1. 



our catalog has 3768 sources detected in J- and Ks-ha,nd 



Number of missing suhstellar objects 

Based on our optical and J K / Spitzer surveys we can 
put limits on the number of very low mass sources that 
are missing in the current census of YSOs in Cha-I. In 
FigureJTJwe show the spatial coverage of our IzJK survey 
(solid lines), the positions of the known members (black 
dots), including the known brown dwarfs (diamonds). 
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TABLE 3 

Candidate members in Cha-I from J K / Spitzer photometry 



ID 


a(J2000) 


<5(J2000) 


J (mag) 


J-K (mag) 


3.6 fivci 


4.5 lira 


5.8 /im 


8.0 Aim 


Notes 


IRl 


11 08 45.38 


-77 33 03.2 


17.956 


2.106 


15.298 


14.912 


14.877 


13.461 


extended? 


IR2 


11 11 23.87 


-76 42 29.1 


18.195 


2.142 


14.847 


14.319 


13.945 


12.916 


extended? 



[Cambresy 
Cnamaeleon- 



the two candidates from Table [s] (squares) , and the two 
new candidate members from spectroscopy (circles) . Our 
fields cover 88 out of 227 known members (38%), from 
which 67 are recovered in our Iz catalog. 64 members 
are recovered in the JK catalog, including 21 with disks. 
Since the object density across Cha-I is highly variable, 
we consider it best to scale the numbers with the fraction 
of known objects covered by the survey, and not with the 
fraction of total area. The scaling factor for obtaining 
the number of potential new candidates in the regions 
n ot covered b y our s urvey is then 1/(1 — 0.38) = 1.6. 

( |1999 ) derived extinction map of the 
-Musca region based on optical star counts. 
Here we define the extent of Cha-I as the region on the 
plane of the sky outlined by a circle of radius r = 1.5° 

and centered at (a - 11 : 05 : 00, 5 77 : 24 : 00), and 

within the contours of Ay > 0.5. Only 1% of this area 
is at Ay > 5. Therefore, when deriving mass equivalent 
to certain brightness limit, we consider only extinctions 
below Ay — 5. 

Let us consider the sources in the candidate selection 
box of the optical part of the survey in Figure[3] that 
are brighter than the limit for spectroscopy, found at 
I « 21. The upper limit placed at / = 17.5 corresponds 
to masses of (0.015 — 0.05) Mq, and the lower limit / = 
21 to (0.004 - 0.008) Mq, in both cases for Ay = - 
5. The estimates are based on DUSTYOO and COND03 
evolutionary models. From the 45 sources in the box, 
15 are known members according to the Luhnian census. 
We obtained spectra for 9 of the remaining 30 objects, 
and find 1 potential new M-type member, falling above 
the substellar temperature limit. From this we conclude 
that there are potentially 2 additional missing members 
in the / = 17.5 — 21 range, in the surveyed part of the 
cluster. The total number of the missing VLM members 
with / < 21 in Cha-I is thus 3xl.6-t-2 = 7 (unsurveyed 
+ surveyed parts of the cluster). 

In our JK/ Spitzer survey, we find zero new objects 
with disks brighter than J = 17 mag, and 2 potential new 
members with disks fainter than J = 17 mag. Accord- 
ing to the COND03 and DUSTYOO models, the J = 17 
hmit corresponds to (0.005-0.008) M© for Ay = - 5, 
similar to the 1 — 21 limit. These numbers have to be 
scaled with the disk fraction (^ 50%, according to Luh- 
man et ar]|2008 1 and the coverage. The total number of 



potentially missing very low mass members in Cha-I is 
< 3 for J < 17 and < 7 for J > 17 mag. This estimate 
relies on three assumptions: 
a) The disk fr action does not chan ge as a function of 
In Cha-I, 



mass 
fraction is rough 



Luhman et al. 
y constant at 



( 2008 ) find that the disk 
SO^rom 0.01 to O.SMp- 



This is not necessarily the case in other star f o rming 



regi ons, as shown, for example, by Lada et al. (20061 



and [Scholz et al.| ( |2007[ ). Below the deuterium- burnmg 
limits, the traction of the disk-bearing objects is even 
more uncertain. 



b) The spatial distribution does not change as a func- 
tion of mass. It is conceivable that brown dwarfs are scat- 
tered towards the outskirts of Cha-I whereas the higher 
mass stars are more concentrated towards the cloud core. 
From the current census, however, there is no evidence 
for brown dwarfs being distributed different than stars 
( |Luhman|[2007) . 

c) i'he Spitzer survey fully covers our JK area. In Fig- 
ure [T] we overplot the area on the deep and shallow IRAC 
images with dashed and dash-dotted lines, respectively. 
Our fields 1 and 4 with the highest density of known 
members (55/88) are completely covered by the deep 
IRAC images. The shallow IRAC images cover about 
5/8 of the remaining fields. However, in these regions 
the density of known members is significantly lower and 
they are offset from the two main cloud cores in Cha- 
I. Therefore we do not consider the missing coverage a 
serious limitation. 

In summary, from both the optical and the 
JK/ Spitzer survey we conclude that the current cen- 
sus of members in Cha-I is mostly complete down to 
^ O.OO8M0 for Ay < 5. Within these limits and barring 
a significant mass dependence in disk fraction or spatial 
distribution, the number of missing object is < 7. There 
may still be a significant number of unidentified members 
below this threshold or at higher extinctions. 

4.2. Ratio of stellar to substellar members 

Based on the current census in Cha-I, we can estimate 
the ratio of low-mass stars (O.OSM© < m < IMq) to 
substellar objects. As substellar objects, we assume all 
objects with Toff < 3000 K, i.e. of the spectral type M6 or 
later. There are 183 low-mass stars in the Luhman cen- 
sus. Here we do not include the two new low-mass-star 
candidates found in this work, since their membership yet 
has to be confirmed. The number of substellar objects 
in the current census is 40. 10 out of these 40 objects 
are of spectral type M9 or later (Tcs < 2400 K). IMyr 
old objects at these low temperatures would have masses 
< 0.015AfQ, according to the DUSTYOO and COND03 
models. Ratio of low-mass-stars to BDs in Cha-I is there- 
fore 183/30 ~ 6, in agreement wi th the 3.3-8.5 limits 



found in other star forming regions ( Andersen et al.|2008 
calculated for the mass range at and above O.OijM©) 
This number represents the upper limit on the star-to 
BD ratio, since we are probably still missing some ob 
jects. 

The ratio between the number of PMOs to BDs is cur- 
rently 10/30 = 33%. This number, however, is subject to 
a large uncertainty, due to missing BDs and also a much 
larger number of missing PMOs. From our statistics, we 
can give a crude estimate of the number of missing PMOs 
in Cha-I. The BD/PMO hmit at 0.015Mq corresponds 
to / w 17.1 - 19.5 for Ay =0-5. For the objects with 
/ < 19.5 in the selection box in Figur e [3} and following 
the same reasoning as in Section |4.1[ we estimate that 
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the number of missing objects is ^ 3. This means that 
the total number of missing BDs is < 3. All remain- 
ing missing objects should be PMOs. Down to / = 21 
[J ^ 17) this gives 4-7 PMOs. From our JK/Spitzer 
survey in Section [4.1| we estimated the number of missing 
objects below this limit to be < 7. In total, this estimate 
gives 4-13 missing PMOs in Cha-I. The PMO-to-BD ra- 
tio in Cha-I then lies in the range between 42% and 77% 
while the ratio of PMOs to the total cluster population 
becomes 5 — 8%. We note that the ratios derived here 
should be taken with caution because they are subject to 
large uncertainties due to scaling from small numbers. 

We can compare the derived numbers with those found 
in oth er star forming regions. In a Ori, Caballero et al. 
( 2007 1 find the lower limit for the PMO versus BD ratio 
of 25 ±1 %. In the ONC, [Lucas et alT] (2006) estimate 
the fraction of PMOs in the total cluster population to 
be in range (1-14)%. Highly variable extinction and age 
uncertainty make an unbiased census in the ONC diffi- 
cult. While these numbers are in line with our Cha-I 
estimates, significantly d ifferent behavior is observed in 



dicators include spatial distribution of stars and BDs in 
star forming regions, possible kinematic differences be- 
tween the two populations, and binary properties (fre- 
quency, separations, mass ratios) and how they depend 
on th e primary mass ( Ahmic et al.|2007 Lafreniere et al.| 
^008|. 



NGC 1333. Scholz et al. (2009) report an overabundance 
of BDs relative to low-mass stars (ratio of 1.5 ±0.3), and 
a deficit of PMOs. They find no PMOs in NGC 1333, 
although clearly being able to do so, as judged from the 
completeness limit of the survey. While the reason for 
this observed difference is not obviouos, it might indi- 
cate that the environment may leave an imprint on the 
shape of the IMF at the lowest masses. The major con- 
tributions to the uncertainties in the PMO /BD ratio are 
introduced by the completeness limits (most of the sur- 
veys are incomplete in the PMO regime), and the mass 
estimates that are based on the models, together with the 
age and distance uncertainties. High and variable extinc- 
tion that is typical for most star forming regions addi- 
tionally complicates this situation, as color-magnitude 
criteria used for candidate selection may also include 
more massive members at higher extinctions. Here we 
must stress the importance of obtaining follow-up spec- 
troscopy of the photometrically-selected candidates, be- 
cause only the spectra can give us the information about 
the youth (and thus membership) and temperature of the 
observed objects. 

As pointed out by [LuhmaJ (2007") , and also apparent 
from Figure [Tl the spatial distribution of brown dwarfs 
does not significantly differ from the distribution of stars 
in Cha-I. This indicates that the ejection from multiple 
systems cannot be the dominant process of brown dwarf 
formation in this cluster. Any other statement concern- 
ing the relative importance of BD formation scenarios is 
not possible based only on ratios of BDs to stars or PMOs 
to BDs. To some extent, any of the proposed models can 
successfully reproduce the observed IMF. For example, 
in simulations of the turb ulent cloud fragmentation by 
Padoan & Nordlund ( 2004 ), abundances of brown dwarfs 
in different star forming regions can be matched by an ap- 
propriate combination of cloud density, temperature and 
the Mach number of the turbulence. To constrain the rel- 
ative contributions of different formation mechanisms to 
the observed substellar population in different star form- 
ing regions, we need secondary indicators that provide a 
more effective ways to distinguish between diflFerent sce- 
narios (see e.g. Bonnell et al.|2007). These secondary in- 



5. SUMMARY AND CONCLUSIONS 

This paper is the third in the series of papers pre- 
sented in the framework of the SONYC project (short 
for Substellar Objects in Nearby Young Clusters). We 
have obtained deep optical and near-infrared observa- 
tions of 4 fields in the Chamaeleon-I star forming region, 
using VIMOS at the ESO VLT and SOFI at the ESO 
NTT. Follow-up spectroscopy of the candidates selected 
from the photometry was obtained using VIMOS/VLT. 
In the following, we summarize the survey and its most 
important results. 

1. In our imaging survey, we reach completeness lim- 
its of 23.0 in /, 18.3 in J, and 16.7 in Ks- This 
corresponds to mass limits of (0.003 — O.OO5)M0, 
for Ay =0—5 and the distance of ~ 160 pc (based 
on DUSTYOO and COND03 models). This is more 
than 2 magnitudes deeper in all three bands than 
the completeness limits of the major part of the 
currently most extensive Cha-I su rvey pr esented in 



Luhman| (|2007[ ). Only one field in |Luhman^ (2007^ , 
covermg 0.22" x 0.28° (equivalent to a quarter of 
our coverage) , did achieve photometric depths com- 
parable to those presented here. 



2. The spatial coverage of our survey is 0.25 deg^. It 
includes 38% of the previously known members of 
Cha-I from the Luhman census. 

3. From the optical photometry, we select 142 objects 
with I — z colors as expected for substellar mem- 
bers of the cluster. 45 of these objects are above the 
/ = 21 limit observed in our spectroscopic follow- 
up, and for 18 of these sources we have obtained 
optical spectra. We have also obtained spectra for 
211 randomly selected objects outside the selection 
box, including several known members. We find 
two objects with spectra consistent with the spec- 
tral type M, but do not have enough information 
to clarify membership. None of the two new po- 
tential members is substellar, as the temperature 
estimate resulting from the spectral fitting are well 
above 3000 K. 

4. Our JK photometric survey combined with the 
Spitzer data yields 2 new candidate members. The 
membership for these objects has to yet be con- 
firmed by spectroscopy. 

5. Based on the results of our survey, while scaling 
for the coverage and disk fraction, we estimate that 
the number of the missing low-mass members down 
to 0.008 Mq and Ay < 5 in Cha-I is < 7, i.e. 
< 3% of the total number of members according to 
the current census. We might, however, still miss 
objects with lower masses, and objects at higher 
extinctions. 
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